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Abstract 



The possibility of charge density oscillations in a finite-thickness two- 
dimensional system is investigated for strong magnetic fields and integer filling 
factors. Using an effective action formalism, it is shown that an oscillatory 
charge density (OCD) is generated in a self-consistent way and is favored en- 
ergetically over homogeneous distributions. It is smooth on the scale of the 
sample thickness and of the magnetic length. The modulus of its wave vector 
is shown to be experimentally accessible. The Hall voltage and the current 
density are shown to oscillate with the same wave vector when a weak current 
is applied. The stability of the charge oscillations against impurity potentials 
is discussed. 
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Two-dimensional electron gas systems (2DEG) subject to a magnetic field normal to the 
plane of confinement exhibit a rich phase structure. This includes the manifestation of the 
integer and fractional quantum Hall effect |1[ at low temperature and the expected Wigner 
crystallization at low densities |2[. 

In this letter, we study a system of interacting nonrelativistic charged fermions in 3+1 
dimensions embedded in a uniform neutralizing background of opposite charge. This system 
is subjected to a uniform perpendicular magnetic field and the particles interact with each 
other through the repulsive Coulomb interaction. The goal of the study is to analyze the 
low energy effective theory describing this system and to explicitly show the possibility of 
charge density oscillations. 

From the study of strictly two-dimensional systems in the absence of magnetic fields 
we know that the ground state is a fluid of constant density 0. Furthermore, it has been 
realized that a magnetic field driven charge density redistribution in the bulk is possible 
as long as no Landau level is partially occupied However, in these studies, the self- 

consistent equations used to describe the charge density do not allow for oscillatory type 
solutions. Here, we show that these are permissible and favored energetically over the 
uniform distribution for real systems with finite-thickness when only the lowest subband is 
occupied. 

The approach we use relies on the effective action formalism combined with a derivative 
expansion technique |||7|]. Here, the effective action is simply defined as the sum of all 
connected one-particle irreducible vacuum diagrams in the presence of any background field 
||. This method has been extremely useful, among other things, in showing the origin of 
the Chern-Simons term in 2+1 dimensions || and in analyzing the phase diagram of planar 
superconductivity jHJ. Here, after removing the four- Fermi Coulomb term in favour of a 



Yukawa-type coupling to an auxiliary field <f), the effective action is obtained by integrating 
out the fermions. The background fields in the expansion take account of the magnetic 
field which gives rise to Landau levels as underlying single particle states and also of a self- 
consistent electrostatic potential. The result of the latter is to shift the center of the Landau 



wavef unctions and thereby make the electron density position dependent. In this way it is 
possible to maintain an integer filling factor throughout the sample despite charge redistri- 
bution. Most importantly, we show explicitly that an oscillatory charge density (OCD) can 
be generated in a self-consistent way as long as each Landau level is fully occupied. It is 
smooth on the scale of the magnetic length £b and of the thickness d. The modulus q of its 
wave vector, in the plane of the 2DEG, is given by 

q~q w =- (1) 

v e 1 

Here q£s <C 1, qd w <C 1, u c = \e\B/m* is the cyclotron frequency, m* the effective mass, 
and eo the static dielectric constant which is assumed spatially homogeneous. The Hall 
potential and the current density both oscillate with the same wave vector. A possible 
consequence of this oscillatory charge density is the destruction of incompressibility in line 
with experimental and theoretical studies of finite-thickness effects in a different context 

The action S, at finite temperature T, describing the interacting finite-thickness 2DEG 
is given by 

S = f dr I d 3 x^ {d T + n-H )^-^- f dr J [^V(x, r) - p] [V>V(y, t) - p 

(2) 



o J ' 2eo Jo 



the first term is the kinetic energy with a parabolic dispersion relation for the fermions; 
H = (-ih\7 + eA) 2 /2m* + U(z), the vector potential is A = (-By, 0,0), (5 = 1/T (k B = 1) 
and /x is the chemical potential. The potential U(z) describes the confinement of the 2DEG 
in the z direction. It can be that of a square, triangular or parabolic well. Finally, the last 
term describes the Coulomb interaction between the charged fermions. 

In order to make S bilinear in terms of the fermion fields, we add to it a term of the 
form 

6S = ~ [dr />x [V0(x, r)] 2 (3) 

07T J J 



and perform the shift 0(x, r) — > 0(x, r) + (ie/e ) J(d 3 y/\x. — y\)ip^(y : r), which cancels 
the Coulomb interaction term in (2). The partition function is obtained by integrating over 
the fermion field ip as well as over the scalar field . Clearly, this procedure has no effect 
since the integration over gives a physically irrelevant constant; however, it allows us to 
cancel the quartic term and thereby make S bilinear in the fermion fields. The linearized 
action reads 

S = J drj (fx |V (d T + n-H - *e0(x, r)) V - ^ [V0(x, r)] 2 + *ep0(x, r)| (4) 

It is important to emphasize that the scalar field <fi in Eq. (4) couples to the charge den- 
sity ip^ip; its fluctuations in the presence of the uniform magnetic field will describe the 
magnetoplasmon collective excitations. 

The advantage of the action (4) is that its bilinear nature in the fermion fields allows 
one to eliminate them exactly. The result of this operation is the normalized effective action 
S cS given by 

S cS = J drd 3 x ^ [V0(x, r)] 2 - *ep0(x, r) J 

-Tr ln(d T + /I - H - ie0)(<9 T + /i - H - V)-\ (5) 

where Tr(- ■ •) denotes the functional trace operation. The potential V(x) in Eq. (5) will be 
specified later. 

In this formalism, a uniform charge distribution corresponds to taking 0(x, r) purely 
fluctuating. Here, we assume a nontrivial vacuum structure corresponding to some average 
electric potential 0o( x > r ) = 0o( x )- This can be obtained by minimizing the effective action 
5S c g/5(j)(x,r) = 0; the result after the substitution e0 o (x) = — *^(x), is the usual Poisson 
equation 



W(x) = 



"- (x ' T| a r + M-go-y |x ' T> 



47TP 2 

■—6p(x). (6) 



Furthermore, the fluctuations around this configuration are described by 
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Scff = / rf 3 x [W(x)f + [ rfrrf 3 x [V0(x, r)] 2 

(7) 



— Tr In 



1 -a T + / x-Ho-V(x)^ (x ' r) 



:n(x,y;u;) 



(9) 



The last term in S c s can be evaluated by a derivative expansion technique 0. In a dia- 
grammatic language, it corresponds to a sum of an infinite number of loops with external 
legs representing the field <fi. The quadratic term in this expansion, defines an induced prop- 
agator for the field 0, and the higher order terms represent self-interaction. Here we will 
restrict ourselves only to the quadratic term giving 

2 

- |s E / <P*d*y<f>(y, -w)n(x, y; u)<p(*, u), (8) 

where II(x, y;a>) is the polarization tensor, represented by the diagram in Fig. (1) and 
uj = 2im/f3 the Matsubara frequency. The ground state energy in the presence of the 
nontrivial electrostatic potential is obtained by integrating out the fluctuations <fi. In terms 
of the potential V(x) and the polarization tensor Il(x, y; uj), the result is 

E G = --^ / rf 3 x [W(x)] 2 + lim i-Tr In [l + 

The evaluation of II(x, y; u) requires the knowledge of the single-particle states of the 
Schrodinger operator H + V(x). Before we proceed into that we exploit the fact that when 
a small electric field is applied in the x-direction, a Hall voltage develops in the y-direction. 
As a consequence, the electrostatic potential that enters the Schrodinger equation depends 
only on the y coordinate. Then the eigenfunctions of the Schrodinger operator take the form 
exp(ik x x)ifje(y — yo)X (z)/ y/L^ with yo = £ 2 B k x , Xq(z) the lowest subband wavefunction in 
the ^-direction and ipe(y) satisfying 

/ h 2 d 2 1 \ 

\2^df 2 + 2 mWcy2 + V{yo + V) ~ £e ' yo ) My) = ° ; (10) 
here, £ is the Landau-level index, and ei >yo the eigenvalue measured relative to the lowest 
subband energy. Provided that V(y) is a slowly varying function on the scale of £b, i- e., 
V (yo) <C muo 2 , we can expand V(y+y ) in a Taylor series and replace it by V(y ) +yV'(y ). 
Then ip E (y - y ) = $ e (y - y' ) and 
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where y' = yo — V ' (yo)/m*u^, and is a harmonic oscillator function. The last two 

terms of Eq. (11) can also be obtained by considering V(y) as a perturbation and using 
non-degenerate perturbation theory. 

As in Ref. [|J we assume that no Landau level is partially occupied and that the initial 
charge distribution is homogeneous in the 2D plane when V(y) = 0. Then the charge density 
Sp(y, z) entering Poisson's equation is given by 

sp(y, ^fE [^(y - y'o) - *?(v - yo)} x 2 (z) = ^-v'\ y )xl{z\ (12) 

- x Lk x 



where v is the filling factor. The last equality follows from a Taylor expansion and the 
assumption 5yo = 2k1 2 b /L x <C £b', it can also be obtained if <&e(y — y'o) is evaluated by 
perturbation theory. In either case the entire Hilbert space is used to obtain Eq. (12) 
as the perturbation V(y) ~ V(yo) + yV'(yo) connects a Landau level with the higher and 
lower levels, e.g., $0(2/ — y'o) * s n °t orthogonal to $i(y — yo). This is in sharp contrast 
with usual treatments, e.g., those of Refs. |l3|,[l4l , which are restricted to the subspace of 
the lowest Landau level wavef unctions. Eq. (11) and the last term in Eq. (12) cannot be 
obtained if we restrict ourselves to this subspace. Notice that Eq. (12) can be written as 
Sp(y,z) = 5p{y)Xl(z) with 5p(y) = evV" (y)/2ntiLu c . 

Alternatively, this redistribution of the electron density will create a potential V(y, z) 
given by 

V(y, z) = --f dy'dz'5p( y ', z>) In [(„ - y'f + (z - z'f] , (13) 
where use has been made of the charge-neutrality condition 



/Ly/2 py r -, 

^**^) = Sfe^(V2)-V(-V2)]=0. (M) 



Notice that Eq. (13) merely expresses the potential produced at a point (y,z) by a strip 
parallel to the at point (y ,z ), with charge density 5p(y ,z). 
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We assume, in line with most experiments, that only the lowest subband in the z- 
direction is occupied. The electrons can then be thought of as extended rod-like charges 
in the ^-direction which are allowed to move only in the xy plane. For these reasons, the 
potential that enters the Schrodinger equation is approximated by the average V(y) ~< 
X (z)\V(y, z)\Xq(z) >. Notice that when the thickness of the 2DEG is zero , the integral 
over z is easily carried out and yields Eq. (12) of Ref. {§]. 

Differentiating Eq. (13) twice with respect to y and using Eq. (12), gives the following 
integral equation 

Spiv) = --^7 / Xl{z)dz [ X%{z>)dz' I"" 12 dyMy') } Z '~f~^~ y L - (15) 



e nhu c J uv J uv J-Ly/2 rv " ' \{z' - zf + (y' - yf] 

An oscillatory charge density Sp(y) = 8p(q)sm(qy), with q > is a solution. When substi 
tuted on both sides of Eq. (15), it gives the dispersion relation that q has to satisfy 
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J dzX*(z) - || Xl(z)dz J Xl(z')dz' e- q \ z '- z \ (16) 



where rj = eohuj c /2e 2 v. An oscillatory charge density of the form Sp(y) = 5p(q) cos(qy) is 
another solution to this integral equation; however, it is not physically permissible since it 
violates the charge neutrality condition (14) which allows only odd solutions. 



We evaluate this dispersion for a square well Xq(z) = J2/dcos(7rz/d), —d/2 < z < d/2. 
With s = qd, we obtain 

2ff2 167r4 

S 2 + 47T 2 S(S 2 +47T 2 ) 2 

When qd <C 1, this gives Eq. (1) with q w = 3/d. For a triangular well we use X (z) = 
(b 3 /2 /V2)zexp(-b z/2) with < z, z < oo. For the average thickness do = 3/6o we obtain, 
for qdo 1, again Eq. (1) with q w = 9/8d . Eq. (1) is also valid for a parabolic well of 
frequency u z , for which X (z) = (l/ir£ 2 z ) 1/4 exp(-z 2 /2£ 2 Z ), with q w = w27r/4- 



This solution for Sp(y) and Eq. (12) give V(y) = —(27ihu c /ei , q 2 )5p(q) sm(qy). For com- 
pleteness, we also give the Hall current density, averaged along the ^-direction. By noting 
that x = u c (y — y ), it is given by 



My) 



euj r 



J2(y - yo)^(y - y'o) 



ev 
2ttK 



V'(y) 



(18) 



where the same approximations have been made as in the derivation of Eq. (12) i.e. V <C hu c 
and [V (yo)} 2 / mu 2 <C %oj c . As was noted in Ref. ||], Eq. (18) merely expresses the fact 
that the current density is such that the Lorentz force is locally canceled by the electrostatic 
potential. Furthermore, by integrating over y we obtain that the Hall current is proportional 
to the potential gradient, and the constant of proportionality is e 2 v/2-K%. 

Now that the electrostatic potential profile is known we are able to evaluate the ground 
state energy from Eq. (9). The first term in this equation is the Hartree energy and is easily 
evaluated. The result, per particle, is 

En TT 2 h 2 



W 



-G 



v 2 q 2 m 



(19) 



where Sn(q) = 6p(q)/e and N e = vL x L y /2ir£ B . The minus sign shows that the OCD state 
is energetically favorable. This conclusion remains valid even when we include the next term 
in the exact expression Eq. (9), here we give only the Fock exchange energy contribution. 
This is obtained by keeping only the first term when we expand the logarithm and pass to 
the zero temperature limit which amounts to replacing (1//3) XL by (l/27r) / duj. The result 



is 



Ccxc 
C G 



N e e N ( 
which after some algebra becomes 
e 2 L 



d 3 Q 1 

(2tt) 3 Q2 



doj Q : n(Q,o;), 



(20) 



Cexc 



2ne N e 



E Jdk x jdk' x J J 



ei'=o 



dQ y dQ z F(Q z )G e ,e>(u) 
(2tt)2 ( kx - K) 2 + Ql + Q 2 z 



(21) 



Here 



u 



Q 2 " 2 



y<-B 



+ f 



G v (u) = (£\ /£'\)e- u u e 
F(Q z ) = \(X \e^\X 



k v k x + 



V'(£ 2 B k x )-V'(£ 2 B k' x y 2 



UJr 



u 



(22) 
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and L\ ~ l is an associated Laguerre polynomial for £' > i (for £' < £, interchange i and 
£' in this formula). This Fock exchange energy is a very complicated functional of the 
potential V or equivalently of the charge density amplitude 5p(q). It can be verified after a 
straightforward algebra that the term linear in Sp(q) is zero. The quadratic term is given, 
for v = 1 and d < £ B 

Sn(q) 



S^ c ( d \ e 2 

G -- 0.615-0.051- ' 



n 



(23) 



N e V ?bJ e £ B 

So far, we have examined the OCD state at zero temperature. It is not hard to extend the 
same analysis to include the temperature effects. In order to do so, one has to evaluate the 
free energy when the charge density oscillates. Denoting by 5F the difference in free energies 
between the normal, spatially homogeneous, and the OCD, spatially inhomogeneous, states 
with fixed N e . This difference can be written in terms of the thermodynamic potential as 
5F = Vt OCD (p) - ^o(Po) + N e (p - p ). Because N e = -dQ /p , 5F = Q OCD (p) - Q (p) + 
l/2(Ap) 2 d 2 Q /dpl H where Ap = p - p . Using V{9) = V sm6, and 

we obtain for low temperatures, and \V \/T ^ 1 

SF = w- T ^ e'^/ T , (25) 
'27H Vol 



where 7 = 1/2 — (Ap + |Vo|)/^o; c . This gives a negative contribution, and hence proves the 
stability of the OCD at low temperature. 

The OCD remains also stable in the presence of weak disorder as is usually the case in a 
high-mobility 2DEG. This conclusion is reached by computing the impurity averaged Green's 
function which satisfies Dyson's equation. Now the relation defining the charge density in 
Eq. (12) is modified by multiplying the squared wavefunctions <&e(y) by the appropriate 
spectral functions. A straightforward calculation reveals that the only effect of disorder is 
Landau level broadening and Eq. (12) is unaffected. 

We now comment on the use of the approximation V(y) ~< X (z)\V(y, z)\X (z) >. It 



has been used in the past to study the collapse of the FQHE state in a quantum well |12[ . It 



is reasonable as long as the electron density is not too high and the quantum well is not too 
wide such that the subbands are not close to each other, otherwise the wavefunction X (z) 
will be modified. Finally, in evaluating the ground state energy in the presence of the OCD 
we limited ourselves to the Hartree-Fock contribution although Eq. (9) goes beyond this 
approximation. This justified by the fact that the neglected terms come in higher powers of 
(e 2 /eo^s)[5n(g)/n] and are small for strong magnetic fields under the assumptions used to 
derive Eq. (12). 

An estimate of the parameters that enter Eq. (1) shows that the OCD is experimentally 
accessible. For instance, for a square well based on GaAS with d = 140A, m*/m = 0.07, 
e = 12 and B — 15 Tesla we obtain A = 2ir/q = 2/j,m. A consequence of the OCD is that 
the Landau levels given by Eq. (11) oscillate with y Fig. I. This facilitates interlevel 
transitions through the indirect gaps and may destroy the quantum Hall effect. 

In summary we have shown that for fully occupied Landau levels and strong magnetic 
fields, a self-consistent oscillatory charge density (OCD) is possible in a finite-thickness two- 
dimensional system, when the lowest subband is occupied. The OCD, which is impossible for 
the ideal two-dimensional system, is smooth on the scale of the magnetic length and of the 
thickness. It is energetically favorable over the homogeneous distribution at low temperature 
and remains stable against weak impurity scattering. 
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FIGURES 

FIG. 1. Polarization diagram (bubble) entering Eq. (8). Solid lines are fermion propagators, 
while dashed lines are the scalar fields. 

FIG. 2. The first lowest Landau levels (£ = 0, i = 1) as a function of yo/^B, obtained from 
Eq. (11) when V(y) = Vos'm(qy) and q£s = 1/3. The dashed and solid curves are for Vq = 0.1Hlo c 
and Vq = 0.45?uj c . The dotted curve shows the Fermi level \i. 



12 



